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ABSTRACT 


The  results  of  an  experimental  investigation  of  the  effects  of  bit 
synchronization  error  and  carrier  frequency  uncertainty  on  the  perfor¬ 
mance  uf  receivers  used  to  demodulate  differential  phase  shift  keyed 
(DPSK)  signals  are  presented.  Also,  previous  analytical  results 
relating  to  the  differential  detector  (DD)  are  summarized  and  an 
expression  for  the  performance  of  the  perfectly  timed  DD  with  frequency 
uncertainty  is  derived.  The  bit  error  probability  and  the  probability  c 
two  consecutive  errors  are  used  as  performance  measures.  The  exper¬ 
imental  and  analytical  results  obtained  apply  to  cases  in  vhich  the 
carrier  is  biphase  modulated  with  the  differentially  encoded  bit  stream, 
the  (equivalent)  noise  at  the  receiver  input  is  white  and  Gaussian  dis¬ 
tributed,  the  bit  synchronization  error  ’S  Gaussian  distributed,  and  the 
error  in  the  carrier  frequency  estimate  is  constant  as  a  function  of 
time.  The  results  presented  provide  a  basis  for  the  design  of  systems 
employing  DPSK  and  differential  detection  when  changing  system  geo¬ 
metries  must  be  considered,  e.g.,  in  satellite  communications  systems. 
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CHAPTER  I 
INTRODUCTION 


The  combination  of  aifferential  phase  shift  keying  (DPSK)  and 
differential  detection  is  known  to  provide  optimum  performance  in  digi¬ 
tal  communications  systems  which  emp’oy  a  carrier  signal  when  an 
a  priori  estimate  of  the  carrier  phase  is  not  available  at  the  receiver. 
However,  to  date,  DPSK  (and  differential  detection)  has  not  teen  utili¬ 
zed  extensively  in  conventional  space  communications  and  telemetry 
systems  because  an  accurate  estimate  of  the  carrier  phase  can  usually 
be  established  and  continuously  maintained  at  the  receiver.  As  a 
result,  phase  shift  keying  tPSK)  anJ  (partially)  coherent  Infection 
usually  provide  best  performance  .  Recently,  considerable  interest 
has  developed  in  non -conventional  systems  in  which  pulsed-envelope 
signals  are  employed  to  c<  nvey  digital  information,  e.g-i  time  division 
multiple  access  (TDMA)  space  communications  systems[J3]  .  In  such 
systems,  it  is  often  either  imo  ssible  or  unpractical  to  retain  phase 
information  accurately  from  pulse  to  pulse.  On  assuming  that  phase 
information  would  have  to  be  acquired  at  the  beginning  of  each  pulse, 
that  the  bit  timing  information  required  to  identify  each  signalling  in¬ 
terval  can  be  retained  from  pu’ae  to  pulse,  and  that  fewer  than  approx¬ 
imately  one  thousand  data  bits  are  to  be  conveyed  pec  pulse,  DPSK  and 
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differential  detection  can  provide  better  overall  performance  than  PSX 
and  (partially)  coherent  detection.  Differential  detection  is  markedly 
superior  to  partially  coherent  detection  in  TDMA  systems  when  only  a 
few  tens  of  bits  are  transmitted  per  puise  and  the  assumptions  just 
stated  are  applicable[5]  .  Thus,  there  is  a  need  for  detailed  informa¬ 
tion  regarding  the  design  and  performance  of  practical  differential 
detectors . 

The  purposes  of  this  report  are  t  reefold:  1)  to  summarize  the 
available  analytical  results  relating  to  differential  detection,  2)  to 
extend  available  analytical  results  by  determining  the  effects  of  an 
error  in  the  carrier  frequency  estimate  on  the  performance  of  the  dif¬ 
ferential  detector  (DO),  and  3)  to  present  extensive  experimental  re¬ 
sults  describing  the  performance  of  a  practical,  i.e.,  imperfect,  DD 
when  the  channel  is  Gaussian.  The  experimental  results  verify  the 
accuracy  of  four  theoretical  results:  1)  the  results  relating  the  bit 
error  probability,  Pg,  to  the  bit  energy  to  single-sided  noise  density 
ratio,  E^/no,  when  the  DD  is  ideally  implemented [  6-9] ,  2)  the  results 
obtained  by  Saiz  and  Seltzberg[l  1J  which  show  the  probability  ol  two  con 
secutive  errors  as  a  function  of  E^/Nq  when  the  BD  is  ideally  imple¬ 
mented,  3}  the  resultJ  derived  by  Huf£[l]  relating  Pg  to  E^/Nq  when 
the  DD  is  imperfectly  timed,  and  4)  the  analytical  results  to  be  derived 
which  show  the  effects  of  an  error  in  the  estimate  of  the  carrier 
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frequency  on  Pg  when  the  DD  is  otherwise  ideally  instrumented. 
Experimental  results  are  also  given  which  show  the  combined  effects 
of  bit  timing  error  and  in  error  in  the  carrier  frequency  estimate  on 
pg  and  on  the  probability  of  two  cor.seciitive  errors,  for  which  no 
corresponding  theoretical  results  are  available.  All  results  apply 
when  biphase  modulation  is  employed,  and  when  the  channel  simply 
attenuates  and  delays  the  transmitted  signal,  the  receiver  bandwidth 
is  much  larger  than  the  signalling  rate,  and  the  (equivalent)  noise  at 
the  receiver  input  is  white  and  Gaussian  distributed. 


The  effect  of  frequency  offset  on  differential  detection  has  been  deter 
mined  independently  by  Henryf  10J  .  The  different  approaches  employ 
ed  in  the  two  analyses  lead  to  expressions  for  Pp  which  differ  con¬ 
siderably  in  form.  However,  spot  checks  of  the  numerical  results 
obtained  using  the  two  expressions  indicate  tnat  they  are  equivalent. 
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CHAPTER  II 

A  DESCRIPTION  OF  THE  DIFFERENTIAL  DETECTOR  AND 
A  SUMMARY  OF  AVAILABLE  ANALYTICAL  RESULTS 


A.  Introduction 

In  this  section,  analytical  results  applicable  to  differential 
detection  are  summarized.  This  summary  is  preceded  by  a  descrip¬ 
tion  of  the  signalling  waveform  and  methods  for  differentially  detect¬ 
ing  the  data  carrying  signal  in  the  presence  of  noise.  Notations, 
assumptions,  and  discussions  applicable  to  subsequent  analyses  are 
given. 


B •  System  Description 

In  systems  employing  (antipodal)  DPSK,  the  carrier  signal  is  bi¬ 
phase  modulated  with  a  differentially  encoded  bit  stream.  Either  of 
two  coding  conventions  can  be  used  to  generate  the  modulation  wave¬ 
form  from  the  data  bit  stream.  To  be  specific,  it  is  assumed  that  a 
"0"  in  the  data  bit  stream  is  encoded  as  (0,  0)  or  (1,  1),  i.e.,  the  en¬ 
coded  bit  in  the  current  interval  is  the  same  as  the  encoded  bit  in  the 
previous  interval.  Correspondingly,  a  "1"  in  the  data  bit  stream  is 
encoded  as  (0,  1)  or  (1,  0).  This  encoding  convention  is  illustrated  by 
exai  lple  m  Fig.  1.  Note  that  the  first  signalling  interval  is  devoted  to 
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SIGNALLING  INTERVAL 
Fig.  1.  A  differential  encoding  example. 

establishing  a  reference,  so  that  N  +  1  intervals  are  required  to  trans¬ 
mit  N  information  bits. 

The  transmitted  waveform  is  assumed  expressible  as 
(!)  SQ(t')  =  Ao  cos[wct’  +  iri(t')  +  q0] 

where  i(t')  represents  the  encoded  bit  stream  and  is  constant  at  either 
0  or  1  during  a  given  bit  interval.  The  amplitude  AQ  is  assumed  con¬ 
stant.  At  the  demodulator's  input,  the  message- carrying  umponent  of 
the  received  signal  is  assumed  expressible  in  terms  of  the  receiver 
time  base  as 

(2)  b(t)  =  s/2p“"  cos  [(wc  +  Awc)t  +  it  i(t  -  «)  +  a] 
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where  Pr  represents  the  power  in  the  received  signal,  «  represents  the 
bit  timing  error,  and  Awc/2it  represents  the  doppler  shift  experienced 
by  the  transmitted  signal.  In  general,  Aw0,  «,  and  a  are  time-varying 
quantities.  An  equivalent  channel  noise  voltage,  n(t),  additively  com¬ 
bines  with  S(t)  at  the  receiver  input  to  form  the  signal  Z(t): 

(3)  Z(t)  =  S(t)  +  n(t) 


It  is  assumed  that  n(t)  is  a  sample  function  from  a  zero  mean  white 
Gaussian  process  with  an  autocorrelation  function  given  by 

<4>  R„nV)=^  N0 


and  a  spectral  density 


(5) 


-CO  <  OJ  <  CO 


Addition' ,  the  noise  process  n(t)  is  assumed  to  be  ergodic. 

A  simplified  block  diagram  of  the  differential  detector  (DD)  to 
be  investigated  is  shown  in  Fig.  2  an  alternative  detector  configura¬ 
tion  is  shown  in  Fig.  3  .  Mathematically,  the  ideal  integrate/dump  dif¬ 
ferential  detector  (i/D  DD)  and  the  ideal  matched  filter/delay  line  dif¬ 
ferential  detector  (mF/DL  DD)  are  identical.  Although  not  proven 
rigorously,  there  is  reason  „o  believe  that  results  to  be  presented  for  an 
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Fig.  2.  A  simplified  block  diagram  of  the 
differential  detector  to  be 
investigated. 


Z(t) 


INPUT 


Fig.  3.  A  second  possible  implementation  of 
the  differential  detector. 
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imperfect  1/ D  differential  detector  apply  when  the  same  imperfections 
exist  in  a  MF/DL  DD[10j  .  Additional  causes  of  performance  degrada¬ 
tion  must  be  considered  when  a  MF/DL  DD  is  employed,  e.g.,  inter¬ 
symbol  interference  caused  b/  imperfect  filtering. 

Referring  to  Fig.  2,  the  equivalent  input  signal  Z(t’  is  first 
separated  into  (approximately)  orthogonal  components  by  multiplying  it 
with  two  phase-quadrature  sinusoids  generated  at  the  receiver.  The 
resulting  multiplier  output  waveforms,  £x(t,  and  £y(t),  contain  both 
baseband  and  double  frequency  components.  Ti  ese  waveforms  are 
processed  by  integrate/dump  circuits.  Ideally,  the  integration  inter¬ 
vals  of  length  Tb  seconds  coincide  with  ihe  intervals  over  which  the 
phase  of  the  received  signal  is  constant.  At  the  end  of  each  bit  inter¬ 
val,  the  integrator  outputs  are  reset  to  zero.  In  practice,  the  inte¬ 
grate/dump  commands  are  not  precisely  synchronized  to  the  bit  arri¬ 
val  time  due  to  bit  timing  error.  Additionally,  practical  integrators 
cannot  be  reset  to  zero  instantaneously.  The  time  required  to  reset 
(or  dump)  the  integrator  will  be  designated  by  6:  a  positivt  constant. 

At  the  end  of  the  m-th  integration  interval,  the  integrator  out¬ 
puts  are  given  by 

pmTb 

exm  *  )  =  ®xm  "*■  nxm 

J(m-l)Tb+6 


and 


rwuj 


(7)  * 


where 


(8) 


and 


(9) 


pirTb 

'  Z(t)4>y(t)  dt  =  Sym  +  Hyjyj 

nm-l)Tb+6  7111 


Sx, 


-  C 


mTb 


ym  -  \  S(t)  <j>Xf  y(t)  dt 

l/(.m-l)Tb+6 


nx,ym=  \  ,v  n(t>  «j>x  y(t)  dt 
y  J(m-l)Tb+6  Y 


Us'ng  straightforward  calculations,  it  can  easily  be  shown  that 
SXm  and  sym  are  given  by 


(10) 


(11) 


^Eh  (^"Tb 

S^J,  =  -= —  \  cos [ it i(t  -€  )]  cos(Aw0t  +  a- y)  dt 

lb  J(m-l)Tb+6 


pmTb 

--  _ _  \  cos[iri(t-«  )]  sin(Aw0t  +  o  -  Y)dt 

7  Tb  ^(m-l)Tb+& 


2  It 


u 


when  »>  ■=—  and  Aw0Tb  <  v  ,  the  case  of  interest.  In  these 
lb  ^ 

expressions,  is  the  radian  frequency  of  the  local  oscillator’1'  (LO), 
Eb  is  the  energy  per  bit  of  the  received  signal,  defined  by 


ir.  most  systerrs,  at  least  one  down- converter  is  employed  prior  to 
detection.  Thus  the  "carrier"  frequency  of  the  signal  being  detected, 
or,  equi1,  ..lently,  the  LO  frequency,  is  actually  an  intermediate 
frequency  (IF). 
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(12) 


Eb  =  PrTb 

and  y  represents  the  random  phase  of  the  L.O.  The  radian  frequency 
Au0  represents  the  error  in  the  estimate  of  the  carrier  frequency  (in 
radians  per  second),  and  is  defined  as 


(13)  A«0  =  +  Auc  -  «d  *  2-rr  A1'0 


where  AfQ  denotes  the  frequency  offset.  It  will  be  assumed  in  subse¬ 
quent  discussion  that  Eqs.  (10)  and  (11)  are  exact. 

The  noise  components  nxm  and  ny^  are  given  by 


(14)  n 


xm 


^mTb 

im  -  I 


n(t)  cos(u>dt  +  Y)dt 


(m-l)Tb+6 


(15) 


pmTb 

nym  “  J 

J(m-1 


(m-l)Tb+6 


n(t)  sin(wdt  +  y)dt 


The  noise  component  are  known  to  be  independent,  zero  mean  Gaussian 
processes  having  equal  variances  given  by [3] 


(16) 


Voltages  exm  and  eym  in  Eqs.  (6)  and  (7)  (approximately)  repre- 

A  A 

sent  the  components  ot  4.(1)  in  the  orthogonal  signal  space  (4>x»4>y).  The 


signals  present  at  the  integrator  outputs  at  time  t  -  mT^  can  therefore 
be  interpreted  as  the  components  of  a  vector  Zm,  where 

—  A  A 

(17)  zmsexm^  +  eym^y  * 

Similarly,  the  message-carrying  component  can  be  interpreted  as  de¬ 
fining  a  signal  vector  Sm: 

(18)  S  =  SVTT1  &  +  S  $ 

'  m  xm  T  °ym  Yy 

In  differental  detection,  decisions  are  based  on  the  dot  product 
of  the  vectors  Zm  and  Z(m_i)*  -At  the  end  of  the  m-th  decision  inter¬ 
val,  the  dot  product  V  is  given  by 

(19)  V=Zm.Z(m.1)=|zm||z(m.1)|co^  , 

where  <>  is  defined  as  the  angle  between  Zm  and  Zm_j .  In  terms  of  the 
integrator  outputs  at  times  ti  =  (m-l)Tb  and  t2  =  mTb,  Eq.  (19)  reduces 
to 

(20)  V  =  exm  ex(m  •  1 )  +  eym  ey(m- 1 ) 

The  value  of  V  is  compared  to  a  threshold  voltage  in  the  decision 
circuitry.  If  the  voltage  V  is  greater  than  a  threshold  Vj,  then  the 
output  of  the  threshold  detector  is  a  positive  voltage;  otherwise  the  out¬ 
put  is  zero,  if  the  a  priori  probability  of  a  "0"  in  the  transmitted  bit 
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stream  is  equal  to  one-half,  it  can  be  shown  that  the  performance  of 
the  ideal  differential  detector  is  optimized  by  setting  V<p  equal  to  zero 
[6,  8]  .  For  this  case,  the  decision  boundary  is  independent  of  the  vec¬ 
tor  magnitudes.  Formidable  difficulties  are  encountered  in  attempting 
to  determine  the  optimum  decision  rule  when  the  DD  is  imperfectly 
timed  and/or  when  an  error  exists  in  the  carrier  frequency  estimate. 
However,  even  if  the  optimum  decision  rule  were  known,  the  decision 
rule  applicable  under  ideal  conditions  would  probably  be  employed  to 
circumvent  instrumentation  difficulties.  Thus,  it  is  assumed  that 
receiver  decisions  are  based  on  the  relations 

(21)  exmex(m-l)  +  eymey(nv  1)  <°  decide  "1" 

and 

(22)  exm  ex(m - 1 )  +  eym  ey (m - 1 )  >  0  decide  "0" 

In  the  absence  of  decision  errors,  the  DD  reproduces  the  data 
bit  stream  --  not  the  encoded  bit  stream  --at  its  output.  Since  the 
received  signal  is  noisy,  some  receiver  decisions  will  be  in  error. 

The  probability  that  a  given  decision  is  in  error,  i.e.,  the  bit  error 
probability  Pg,  will  be  used  as  one  measure  of  performance.  The 
probability  that  two  errors  are  made  in  succession,  Pjtj-j,  use<3 

as  a  second  performance  measure-  Analytical  results  summarizing 
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the  performance  of  an  ideal  DD  are  contained  in  the  following  r-' 
section. 


C.  Ideal  Differential  Detection 

1 .  The  bit  error  probability 

The  bit  error  probability  (BEP)  of  the  ideally  intrumented  dif¬ 
ferential  detector  is  given  by  [  6-9] 

(23)  PE  =  j  exp  [  -  Eb/N0] 


where  Eb/No  represents  the  bit  energy  to  noise  density  ratio  of  the 
received  waveform.  As  shown  in  Appendix  A,  the  ideal  DD  performs 
nearlv  as  well  as  the  ideal  coherent  detector  at  high  bit  energy  to  noise 
density  ratios.  When  the  integrator  dump  time  is  non-zero,  then  the 
BEP  of  the  otherwise  ideally  implemented  DD  is  given  by [  1] 


(24) 


PE6  =  j  exp 


Pr(Tb-&) 

No 


2 .  The  probability  of  two  consecutive  errors 

An  analysis  of  the  paired  error  probabilities  of  the  ideal  DD  was 
performed  by  Salz  and  Saltzberg[  1 1]  which  applies  when  the  channel 
noise  is  a  zero  mean  white  Gaussian  process.  First,  an  expression 
for  the  probability  of  two  consecutive  errors,  PEq»  was  found.  This 
result  was  then  divided  by  the  bit  error  probability  to  obtain  an 
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expression  for  the  probability  of  error  given  that  an  error  occurred  in 
the  previous  bit  interval; 


(25)  Fr(emle/^_i))  =  — — 

In  a  form  convenient  for  numerical  evaluation,  the  expression  for 
*Mem !  e(m- 1 ))  as  a  *unct*on  the  bit  energy  to  noise  density  ratio 
was  found  to  be  given  by 

\26)  Pr(em|e(m_l))  =  JL  \  a2(9)[l  r'^IZicos  0exp(Ma  006*0 )b(0)] 

2ir  Ja  ja 


wnere 


a(8)» 


•it 


exp(  -  u2  )  du 


cos  0 


b(0)=  — —  C  exp(-u2)du 

J-Mcos0 


The  detection  errors  associated  with  non- adjacent  symbols  are  inde¬ 
pendent  when  the  DD  is  ideally  instrumented.  Thus,  the  conditional 


probability  Pr(emle(m.n))  for  n  >  2  is  approximately  equal  to  the  bit 
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error  probability  Pg.  As  a  result,  the  occurrence  of  more  than  two 
consecutive  errors  is  a  relatively  rare  event. 

The  results  of  a  numerical  cr  -nputation  of  Eq.  (26)  are  given  in 
Fig.  4.  These  results  show  that  the  ratio  Pr(em  ]  e(m-l))  goes  to  zero 
as  the  bit  energy  to  noise  density  ratio  is  increased  without  bound. 

This  behavior  is  best  explained  by  Masonson(4] :  "An  isolated  single 
error  can  occur  at  very  high  signal  to  noise  if  two  deviations  are  just 
large  enough  to  differ  by  more  than  the  half  angle  p  ,  but  not  so  large 
as  to  imply  that  another  error  will  necessarily  follow.  For  example, 
in  differential  binary  systems,  an  error  occurs  if  two  deviations  suc¬ 
cessively  are  -45°  and  -46°.  In  low  noise,  the  next  phase  deviation  is 
almost  sure  to  be  near  zero  and  so  a  second  error  does  not  occur." 
This  statement  is  further  substantiated  by  the  experimental  results 
given  in  Section  IV-  0. 

D.  The  Imperfectly  Timed  Differential 

Detector 

The  effect  of  bit  timing  error  on  the  bit  error  probability  of  an 
otherwise  ideal  DD  has  been  calculated  by  Huff[l]  .  Explicit  numeri¬ 
cal  results  were  reported  for  the  case  where  the  timing  error,  t,  is 
Gaussian  distributed  with  mean  value  «  and  standard  deviation  <re  .  The 
loss  mechanism  when  the  timing  error  is  non-zero  will  now  be  deicri- 
bed  and  appropriate  results  summarized. 
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Pr<e2|e,) 


o.eo 


Fig.  4.  The  probability  of  error  given  that  an  error  was 
detected  in  the  previous  bit  interva’  versus  the 
bit  energy  to  noise  density  ratio  for  the  ideally 
instrumented  DD. 
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If  it  is  assumed  that  the  frequency  offset  and  the  integrator  dump 
time  are  identically  zero,  and  that  p=a-y  changes  only  a  few  degrees 
during  two  consecutive  bit  intervals,  then  Eqs.  (10)  and  (11)  reduce  to 


(30) 


VC  fTb 


3xl 


Tb 


f  b 

\  cos[’ri(t-«  )]  dt  cos  p 


(31) 


sy*=  “Tb" 


PTb 


c  ° 

\  cos  [iri(t  -  «)]  dt  sin  P 

O  n 


VET  czz b 


_  N/Eb  C 

(32)  Sx2  cos[7ri(t-«)]dt  cosp 


and 


(33) 


y? 


VC  rZT b 

~tT 


1. 


cos[iri(t  -  c  )]  dt  sinP  . 


Correspondingly,  Si  and  are  given  by 

_  -jrb  rTb  a  a 

i34j  Si  =  j  cos  [ni(t  - «  )J  dt  [cos  p<j>x  fsinp4>y] 


and 


Fcr  convenience  of  notation,  m  is  set  equal  to  1  with  the  understand¬ 
ing  that  each  integral  is  performed  over  a  bit  interval  that  is  interior 
to  the  transmitted  bit  stream. 
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(35) 


Sz  = 


Tb 


r 


cos  [iri(t  -  «)]  dt  [cos  P4>x  +  sin(3$y] 


In  the  absence  of  noise,  the  bit  timing  error  affects  the  output  of  the 
integrator  as  shown  in  Fig.  5.  This  figure  illustrates  that  the  sample 
value  Sxm  depends  on  the  bit  sequence,  the  magnitude  of  e,  and  the 
sign  of  «  .  When  a  phase  transition  occurs  in  the  received  wavefoim 
during  the  m-th  integration  interval,  the  bit  timing  error  causes  a 
reduction  in  the  magnitude  of  Sxm .  Correspondingly,  the  signal  to 
noise  ratio  at  the  integrator  output  is  reduced  since  the  statistics  of 
the  noise  components  nx>  ym  are  independent  of  the  timing  error.  When 
no  transition  occurs  in  the  integration  interval,  the  signal  magnitude 
remains  unaffected,  provided  the  bit  timing  error  does  not  exceed  the 
bit  length.  Table  1  shows  the  effect  of  bit  timing  error  on  the  magni¬ 
tudes  of  the  signal  vectors  Sj  and  Sj,  L,j  and  L»,  for  each  of  the  six¬ 
teen  possible  four  bit  sequences  {io,  ij,  i*t  •  As  is  evident  from 
Table  1,  the  probability  that  an  error  is  made  in  detecting  the  infor¬ 
mation  conveyed  by  |  i2  -  ij  depends  on  the  bit  timing  error  and  the 
particular  sequence  {io,  ij,  i2,  13}.  Clearly,  the  probability  of  error 
when  the  data  bit  is  a  "1",  Pj,  is  larger  than  the  probability  of  error 
when  the  data  bit  is  a  "0",  P0|  when  the  timing  error  is  non-zero. 

This  is  a  consequence  of  the  fact  that  the  signal  phase  experiences  a 
transition  between  bits  when  a  "1"  occurs  in  the  data  bit  stream.  The 
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Fig.  5,  An  example  illustrating  the  effects  of  bit  tuning  error  on  the  integrator' 
output  assuming  no  channel  noise  is  present. 


average  bit  error  probability,  Pret,  depends  on  the  relative  frequency 
with  which  the  sixteen  possible  four-bit  sequences  occur.  Huff  has 
shown  that  can  be  expressed  as 

1  n'ffi b* 

(36)  Pe,  ■  -  \  \  K(LlfL2)P(Li.L2)dLidLa 

where 

l2Vnq  d0 

Ll*  cos20  +  L22  sin20_ 

and  P(Li,  L2)  is  the  joint  probability  density  function  for  the  variable 
Lj  and  L2.  Similar  expressions  were  derived  for  the  conditional  bit 
error  probabilities  P0  and  Pt  .  The  expressions  were  numerically 
evaluated  for  the  special  case  where  the  sixteen  possible  four-bit 
sequences  are  equally  probable,  i.e.,  each  sequence  occurs  with  a 
•/robability  of  one -sixteenth.  Huff  assumed  1)  that  the  timing  error 
process  is  independent  of  the  noise  processes  n^^  and  2)  that 

the  tuning  error  is  slowly  varying  with  respect  to  the  data  bit  rate, 
and  2)  that  the  area  under  the  tails  of  the  Gaussian  distribution  function 
oi  the  bit  timing  error  for  |c  |  <  Tt>/2  is  small  compared  to  the  bit 
error  probability,  i.e.,  that  Lj  and  Li  are  essentially  confined  to  the 
interval  (0 , 

Sample  results  showing  P£4,  Po,  and  Pj  as  a  function  of  Et>/N0 
for  «n  »  0  and  <r{n  =  0.08,  where  «n  =  «/T^  and  tr(n  a  (rt/T\3,  are  given 


(37) 


K(Lj,  L2)  = 


pTT/  0 

=  4  exp 


21 


in  Fig.  6.  The  theoretical  results  can  also  be  presented  in  terms  of  a 
normalized  ratio  (3/f3o»  where  po  and  {3  are  defined  as  the  bit  energy  to 
noise  density  ratios  required  to  maintain  a  given  bit  error  probability 
under  conditions  of  perfect  and  imperfect  bit  synchronization,  respec¬ 
tively.  Figure  7  shows  the  values  of  the  ratio  p/0o  as  a  function  of  <r€n 
for  several  different  bit  timing  offsets,  and  for  Prec  *  10~5  .  Note  that 
a  small  bit  timing  jitter  (tr€n  ~  0.01,  «n  =  0)  causes  performance  to 
degrade  less  than  3%  (in  terms  of  (3/Po);  in  most  cases  such  a  loss 
would  be  considered  insignificant.  However,  a  moderate  bit  timing 
jitter  (<r<n  ^  0.06)  causes  performance  to  degrade  more  than  35%  --  a 
relatively  large  loss. 


ZZ 


Eb/N  (  NUMERIC  ) 
”0 


Fig.  6.  The  error  probabilities  Pg,  p0,  and  P 
versus  the  bit  energy  to  noise  density 
ratio  of  the  imperfectly-timed  DD; 
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Fig.  7.  The  ratio  p/,£o  versus  the  normalized  timing 
jitter  for  selected  values  of  the  bit  timing 
offset;  Pj£  =  10"5  . 


CHAPTER  III 

THE  EFFECTS  OF  FREQUENCY  UNCERTAINTY 
ON  THE  BIT  ERROR  PROBABILITY 


A.  Introduction 

In  this  section,  an  expression  for  the  bit  error  probability  as  a 
function  of  the  frequency  offset,  Af0,  and  the  bit  energy  to  noise  den¬ 
sity  ratio,  with  <  =  0  and  6  =  0,  will  be  derived.  Much  of  the  analysis 
parallels  the  approach  used  by  Huff[  1]  in  his  analysis  of  the  imper¬ 
fectly-timed  DD. 


B.  Preliminary  Signal  Analysis 

When  the  phase  angle  a  -  y  [  see  Eqs.  (10)  and  (11)]  changes  only 

a  few  degrees  in  adjacent  bit  intervals,  the  frequency  offset  does  not 

$ 

change  with  time,  and  «  =  6  *  0,  Eqs  .( 10)  and  ( 1 1 )  become 


(38) 


2^oTb 


[sin(Trii  +  Aw0Tb +  P)  -  sinfirii  +  (3)  +  sin(vij-p) 

-  sin(Trii  -  AwQTb  -  P)J 


and 


\lEh 

(39)  Syi  =  -  -  [cos(7rii  + (3)  -  cos(:rii  +  Au>  Tb  +  p) 

2Aw0Tb 

-  cos(irii  -  Aw0Tb  -  (3)  +  cos(irii  -0)] 


Again,  m  is  set  equal  to  1  for  convenience  of  notation. 
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where 


(40)  (3  =  a  -  V  • 


Using  these  expressions  in  conjunction  with  Eq.  (18)  and  employing 
the  appropriate  trigonometric  identities,  the  expression  for  Si  reduces 
to 

(41)  S,(A«0)  =  Li  jcos^P  +  ^“^)#x  .  sin^p  $y 

where 


Lj  = 


sin(AwQTb/ 2) 
A%Tb/2 


“  A  A 

The  relation  of  Si  to  the  -^i^y  coordinate  frame  is  shown  in  Fig.  8. 
The  following  relations  are  readily  derived  from  Eq.  (41): 


Fig.  8.  The  relationship  of  Si  to  the  $x,  $y 

coordinate  frame  with  frequency 
offset . 


(43) 

Iim  Si(AwQ) 

Aw0~-  0 

(44) 

ISi(A«0)|  =  Li 

and 

(45) 

/C(Aco0)  ~/sx(0) 

c.m\  -  AwoTb 


radians  (counterclockwise 

positive) 


A  new  coordinate  frame,  (x',y'),  in  which  Sj  is  aligned  with  the 
x'  axis  can  be  obtained  by  rotating  the  (4>x»  4>y)  frame  (~&-&<x>0Tyy/ 2  +  trij) 
radians  in  the  counterclockwise  direction.  In  ,nis  new  frame  of  refer¬ 
ence,  Si  can  be  expressed  simply  as 

(46)  Si  =  Li  , 

where  3^1  is  a  unit  vector  in  the  plus  x'  direction.  By  employing  the 
same  assumptions  used  to  derive  Sxl  and  SVJ ,  the  components  of  S2  in 
the  original  reference  frame  are  found  to  be  given  by 


(47)  SX2  =  Li  cos  th2  cos(3Aw0Tb/2  +  (3) 

(48)  Sy2  =  -  Li  cos  Triz  sin(3Aw0Tt>/2  +  f3) 
and  thus 

(49)  Sz  =  Licos  n-iz[cos(3A45Tb/2 +p)^x  .  sin(3Aw0xb/2 +|3  )by] 
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In  the  (x*,  y')  coordinate  frame,  S2  is  expressible  as 
(50)  S,  =  ^  cosir  6-  (cosAw0Tb  a x  •  -  sinAw0T|j  ayO 


where 

(51) 


6i  =  Uz  -  il 


l  if  the  data  bit  is  a  "1" 


1^0  if  the  data  bit  is  a  "0" 


Through  inspection  of  Eqs.  (46)  and  (50)  it  can  be  see  that 
(52)  |S, |  =  |S, |  L, 


and 


(53) 


A  ~/s>  =  "  AwoTb 


That  is,  the  frequency  offset  causes  the  magnitude  of  the  signal  vectors 

sin(A«0Tv,/2) 

to  be  reduced  by  the  factor  _ t -  and  introduces  a  Dhase  chut 

hP 

between  and  Si  .  Neithei  the  frequency  offset  nor  the  coordinate 
rotation  affect  the  statistics  of  the  noise  components  and  thus  nx*i, 
nxls  »  ny  and  ny^  are  samples  from  independent,  zero-mean  Gaussian 
processes  having  variance  equal  to  N0/2. 

C.  The  Expression  for  the  Bit  Error  Probability 

In  the  (x1,  y1)  frame,  the  vectors  Zi  and  Z2  are  expressible  as 


—  A  A 

zi  =  ex'i  ax'  +  ey'i  ay 


2  «S 


(54) 


and 

—  A  A 

{55}  Z2  -  ex<i  a*'  +  eyi  ayi  , 

where 

(56)  ex'i  =  Li  +nK'i 

(57)  ey^  =  ny'i 

(58)  ey'2  =  -  Ljcos  tt  5^  sinAw0Tb  +  ny'2 
and 

(59)  exi  =L1coarr6i  cosAo0Tb  +  n^. 

Sample  values  of  Zj  and  Z2  for  the  case  ii  -  i2  are  shown  in  Fig.  9. 

The  angle  between  the  vectors  Sj  and  Zj,  measured  in  the  counter¬ 
clockwise  sense  (see  Fig.  9),  is  defined  as  4* )  •  Similarly,  the  angle 
between  S2  and  Z2  is  defined  as  ij>2 .  For  this  particular  example, 

1 4*i  +  AwoTb  -  4*2 1  <  J  •  thus,  Zi  *  Z2  >  0  and  the  receiver  decides  cor- 

u 

rectly  that  a  zero  was  transmitted.  Since  the  probability  density 
functions  of  4>i  and  iJj2  are  even,  frequency  offset  introduces  an  asym¬ 
metrical  dependence  of  the  bit  errcr  probability  on  the  noise  which 
degrades  performance. 

Assuming  that  the  threshold  voltage,  Vj,  is  zero,  the  error  pro¬ 
babilities  Pq  and  Pi  are  given  by 
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(60)  Po  =  Prob.  i7;  '  z*  <  °/6i  ‘  0) 
>61)  Pi*  Prob.  !Zi-7i>0/6i3  " 


An  expression  for  R  will  now  be  derived. 

Since  n„',  and  nyl  are  statistically  independent  samples  from  aero 
mean  Gaussian  processes,  each  with  variance  N0/2.  the  conditional 
joint  density  function  for  ex',  and  cyl  can  be  written  as  [see  Eqs.  (56) 

and  (57)j 
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(62) 


P  (ex1!*  ey'i  l^/Eb.  Aw0)  =  p  (ex1!,  eyi  |  Li)  = 

JL  exp[~-  ^x*  Jr!1.)2  4e)—  1  . 

irN0  L  N0  J 

The  signal  components  ex*i  and  e^j  can  be  written  in  terms  of  Z*  and 
ij/i  as 


(63)  ex',  =  Uii  COS 

(64)  ey'i  =  |  Zx  |  sinifii 


-  ir  <  i|>i  <  tt 


The  conditional  probability  density  function  of  4* i*  given  Lj,  is  found  by 
substituting  Eqs.  (63)  and  (64)  into  Eq.  (62)  and  then  employing  a  well- 
known  transformation  technique[8]  to  obtain 


(65) 

P(«I»i|Li)  =  -L  exp 

T  2 

+  -4—  exp 

hiZ  sin2ijji 

2tt 

L  N0  J 

\f2rt 

No 

N, 


COS  4*1 


J+  Q(L„«M 


where 


(66) 


Q(L*  i}j) 


_l  r  >iN 

"s/2tt  4) 


COS  ((>1 


exp 


[4] 


du 


Next,  assuming  tpi  given,  a  new  (x,  y)  coordinate  frame  is  formed  by 
rotating  the  (xSy*)  frame  in  the  counter-clockwise  direction  so  that  the 
signal  vector  Zi  is  aligned  with  the  positive  x  axis.  Thus,  in  terms  of 
the  (x,  y)  frame, 


(67) 


eyj  =  (LxCosAw0T|jCosiJj1  +  LisinA«0TbSin4<i)  costtB^  +  ixjj^ 


(68)  ey2  =[-  LxsinAw0Tb  cos  i^t  +  Lx  cos  AwQTb  sini^ijcos  tt  6^  +  ny* 
which,  for  0^  =  0,  reduces  to 

(69)  ex2  =  Lx  cos(Aw0Tb  -  4»x)  +  rixZ 

(70)  ey2  =  Li  sin(Aw0Tb  -  44)  +  ny2 
where 

(71)  nx2  =  nxfe  cos  4>i  +  ny2  sinijJi 

(72)  ny2  ::  ny!j  cosi(ji  -  nx^  simj>i 

As  before,  the  n^y{  are  samples  from  statistically  independent,  zero- 
mean  Gaussian  processes,  each  with  variance  NQ/2. 

Examination  of  the  vector  components  of  Z2  and  Zx  in  the  (x,"y ) 
frame  reveals  that  is  proportional  to  the  dot  product  Zj  •  Z2.  Thus, 
ex2  =  0  corresponds  to  the  decision  boundary  1^2  ■  A<*>0Tb  -  i^i!  =  ^  . 
Therefore,  Eqs,  (60)  and  (61)  reduce  to 

(73)  P0  =  Prob  (exi  <  0|Li,  BiS  0) 

(74)  Pi  =  Problem  >  Oki.  6i  =  i) 


32 


From  Eq.  (67),  the  conditional  density  function  for  can  be 


written  as 


(75)  p(exz|Aw0,  =  0)  =  — - -  exp  -  — [eX2  -  Li  cos(Au>0Tb-iWl  2 

n/ttNq  L  No  J 

By  the  chain  rule  for  marginal  and  joint  density  functions[  12] , 

(76)  P(ex2»  Lj,  4<i  |awo,6^=0)  =  p(e*2|  Li»  A w0,  44,6 i=0) 

*  p(4»J  Li,  Au,0l  6i=0)  p(Lj | Aw0,  6|*0)  . 


Integrating  both  sides  of  Eq.  (76)  over  the  range  of  tjji  and  Lj,  it  is 
found  that 

(77)  p(e-2|Awo,  6i=0)*(  ^  pie**  I  Ll*  A<*>0,  i|q, 6  fO) 

Li  4>l 

•  P(^i|Li,Awq,  6i=0)p(Li|Awo,  5i=0)  di|j1dL1 


For  a  given  Aw0  and  for  6^  =  0,  Li  is  a  constant.  Thus 

1  <=  c!"  1 —  /ain(AwQTb/2)  \ 

(78)  p(Li|Au0,  ©t=0l  =6  L,  -  -JTb  I  A„0Tb/2  '  ) 


where  6(x)  is  the  Dirac-delta  function.  Moreover,  the  value  of  is 
independent  of  the  frequency  offset  and  the  phase  of  the  succeeding  bit 
interval.  Therefore, 


p(«|q  I  I-i.  A«0,  6^0)  =  p(i|q|Lj)  . 
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Utilizing  the  relations  in  Eqs.  (78)  and  (79),  Eq.-  (77)  reduces  to 

(8®)  p(exzl  Lj,  Aw0,  6(=0)  =  \  p(ex2|  Li,  ®i=0)p(^xi  Lj)dij<i  • 

% 

Substituting  Eq.  (80)  into  Eq.  (73),  the  probability  of  error  given  that 
a  "0"  was  transmitted  becomes 


(81) 


f° 

Po  p(exJ  Lj,  A<*>0,  ■<  ;~0)  d  ex2 


-00 
,ir  ^0 


=1.3 


- 1  ~-oo  rN0 


exp 


1  r  _  2 
■  — -  L  ex2  -  L,  cos(Ao>0Tb-ih)J 

wo 


Ide^is 


p(+i  I  Li)#x 


*  J  ’  Q(Ll*AcJoTb  "  •mJ  p(+il  L|)#l 


Similarly,  Pj  is  found  to  be  given  by 


(82) 


r°° 

=  \  p(  6x2 !  Ej,  A  wol  6^3 1)  dejr2 

Jri 


-n 


-tt  0  V* Nq 


exp 


1  r  2 
‘  rpl  eX2  +  Lx  cos(A  w0Tb  -  ipx)J 
1No 

.  p(4»il  E|)d4»i 


dex2 


By  an  appropriate  change  of  variables,  Po  and  Pi  can  be  shown  to  be 
equal;  thus  the  bit  error  probability  is  independent  of  the  transmitted 
symbol.  In  the  remainder  of  the  analysis,  it  will  be  assumed  that 

&i  ■  0. 

Now,  on  substituting  Eq.  (65)  into  Eq.  (81),  <he  expression  for 
the*  total  bit  error  probability  can  be  written  as 
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& 


(83) 


PE  =  P0  -  J-  f  exp  [-  iull  #l  +  f  expLiiWjf] 

4*  J-ir  L  Noj  WiN 'J.n  [  No  J 


cos  i}*«  di|i| 


'  ^  J^Q(Ll*Au,°'rb-^i)exp[' 


2n/ttN0  J.n 


n* 

U  f 

+  2n/ttN^  -v 


f*  r  Lt2  sin2d/i 

)  Q(Li,  Aw0Tb-4<l)  exp  I- - jj- - J  cos  dth 


Q(Lt,4<l)  exp 


£.  L»i2  sin2tjj|  J 


cos  ipj  dip  i 


—  f  Q(Li.4»i)Q(Li,A«0Tb-i(ii)  expf-i^^l 

'VtFnT  J-TT  L  Nq  J 


cos  ip i  dip  i 


Fortunately,  this  equation  can  be  simplified  by  utilizing  symmetry  pro¬ 
perties  of  the  integrands.  It  can  be  shown  that 


(84)  J  Q(L1,Aw0Tb-4'l)  exp  dipi  =  0 


(85)  ^  Q(Li,Aw0Tb-^i)Q(Li»^i)  exp  £- 


r  Li2sin24ii~l  ,  j,  n 

J  exp  -  — — —  J  cos^idi|>i=  0 


L|2  sin2  4^i 
Nn 


cos  i|<i  dipi  =  0 


(86) 


and 
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(87) 


J  Q{Lj,  4<i)exp  -Ll  Jin  cos  4i»  di^i 


pir/2 

=  4  \ 

J0 


Q(L>1*  4u)exp 


J"  L|2sir.2^|1 

L  ng  J 


cos  ipi  dij.'i 


The  expression  for  Pjr  in  Eq-  (83)  is  now  expressible  in  the  form 


(88) 


PE=i  +  E''E‘ 


where 


(89) 


and 


(90) 


E.  = 


2L 


pu/2 

ZZ  \  Q(Li,4ji)  exp 

r 


sltifC  J0 


r  t  2  .  ; 

L |  sin 
N0 


iiil 


cos  ij^i  dij^i 


E2a — — *■  C  Q(I-a,A(d0Tb-4'i) exp  8^n  >l<t~l  cos^id^i  , 

Z\fnW~  -it  L  NQ  J 


It  is  easy  to  show  that[l] 


(91) 


El  =  i-i  «p[-^-] 


The  expression  for  E2  is  evaluated  in  Appendix  B  to  give 

yn7 2 


(92) 


E2 


.  1  1  f1  T  UZ  L22/N0  1  ,Q 

“  — - \  exp - f - V—  9-  - 1 —  d6 

2  it  Jq  L  Ei  cos  0+  L22sin  0 


du 


„i7/2  „L»ljpcos(’I+A“oTb)  r  „ 

+  I  ,  I  rr°  exp  r.f-1 

-tt/2  Lj  — cos  x  cos  Aw0T|j  *■  -* 

^No  , 

•  exp  -  — s  — r-- -x  1  cosx  dx 

L  N0  J 
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where 


(93)  Lz=  Li  cosAco0Tb 


The  desired  expression  for  the  bit  error  probability  as  a  function 
of  the  frequency  offset  AfQ  and  the  bit  energy  to  noise  density  ratio 
Et>/N0  is  obtained  by  combining  Eqs.  (88),  (91)  and  (92).  In  a  form 
convenient  for  numerical  computation,  the  result  is 

,  rW2  r  L,2L22/Nq  "1 

(94)  PE(L„  Eb/N0.  Aw0)  =  -  ^  exp  |-  d® 


•  exp  j-  L‘*ain‘x  j  cosx  dx  . 

The  value  of  the  first  integral  in  the  above  expression  depends  on 
the  projection  of  Sm-1  onto  Sm  and  has  the  same  form  as  the  function 
K(Li,  L2)  defined  by  Huff  [1]  (see  Eq.  (37)).  The  second  integral  takes 
into  account  the  asymmetrical  dependence  of  the  bit  error  probability 
on  the  noise  introduced  by  frequency  offset.  Note  that  the  result  is  in¬ 
variant  to  the  sign  of  the  frequency  offset  and  that,  in  the  limit  as  A«c. 
g.)es  to  zero,  Eq.  (94)  becomes  the  expression  for  the  bit  error  proba¬ 
bility  of  the  ideal  differential  detector,  given  by  Eq.  (23). 


D. 


Ntimerical  Results 


The  expression  for  the  bit  error  probability  was  evaluated 

numerically;  the  results  are  given  in  Fig.  10  for  various  normalized 

frequency  offsets,  Afn  =  Af0Tb»  as  a  function  of  the  bit  energy  to  noise 

density  ratio.  These  results  can  be  presented  in  a  different  manner  by 

plotting  the  additional  signal  power  required  to  maintain  a  given  bit 

error  probability  as  a  function  of  frequency  offset.  Calculated  values 

of  JL  for  selected  bit  error  probabilities  are  shown  in  Fig.  11,  where 
Ho 

pQ  and  p  have  been  defined  as  the  bit  energy  to  noise  density  ratios 
required  to  maintain  the  given  bit  error  probability  for  zero  and  non¬ 
zero  frequency  offsets,  respectively.  The  results  indicate  that  a  small 
frequency  offset  does  not  significantly  degrade  performance,  but  that 
efficiency  drops  sharply  for  larger  offsets.  For  example,  an  offset  of 
(0.01)Tk  results  in  an  effective  energy  loss  of  less  than  0.1  dB  and 
could  be  tolerated  in  most  applications.  An  offset  of  only  (0.05)Tb» 
however,  introduces  a  loss  greater  than  1  dB,  which  may  be  unaccept¬ 
able  . 
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Fig*  10.  The  bit  error  probability  versus  the  bit  energy 
to  noise  density  ratio  for  selected  values  of 
frequency  offset. 
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Fig.  11.  The  ratio  \i/ \i0  (see  p.  38)  versus  frequency 
offset  for  selected  values  of  the  bit  error 
probability. 
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CHAPTER  IV 

THE  EXPERIMENTAL  DIFFERENTIAL  DETECTOR 
AND  ASSOCIATED  TEST  CIRCUITS 


A.  Introduction 

An  experimental  receiver  was  instrumented  to  test  the  validity  of 
the  analytical  results  and  to  provide  additional  results  for  cases  which 
are  extremely  difficult  to  investigate  analytically.  In  what  follows,  the 
experimental  DD  and  the  associated  test  circuitry  are  describe-’,  in 
some  detail. 

B.  Description  of  the  Differential  Detector 

A  block  diagram  of  the  experimental  DD  is  shown  in  Fig.  12- 
The  input  signal  Z(t)  represents  the  received  waveform  at  a  30  MHz  IF. 
This  signal  is  first  split  into  two  components  having  equal  magnitudes 
and  phase  using  a  hybrid  divider.  Two  switch  type  balanced  mixers, 
driven  by  quadrature  local  oscillator  signals,  are  used  to  down-con¬ 
vert  these  components  to  generate  two  (approximately)  orthogonal  base¬ 
band  signals,  £x(t)  and  §y(t).  Each  mixer  is  followed  by  a  low-pass 
gain  stage  which  has  a  bandwidth  equal  to  about  1  MHz.  In  turn,  the 
output  of  each  gain  stage  is  applied  to  an  integrate/dump  "ircuit.  The 
dump  time  of  eaUi  integrator  is  continuously  variable  from  350  nsec  to 
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5  psec.  However,  all  results  to  be  presented  were  obtained  with  6  set 
equal  to  350  nsec,  A  detailed  schematic  diagram  of  the  integrator  cir¬ 
cuitry  is  given  in  Appendix  C. 

The  detector  contains  six  separately  controlled  sample/hold 
(S/H)  modules,  three  in  each  orthogonal  channel.  Each  S/H  module 
samples  the  appropriate  integrator's  output  at  the  end  of  every  third 
interval  and  holds  it  for  two  integration  intervals,  i.e.,  for  2Tb 
seconds.  Thus,  each  S/H  circuit  tracks  for  Tb  seconds  and  holds  for 
2Tb  seconds.  The  three  S/H  modules  in  each  channel  allows  each 
integrator  output  to  be  sampled  at  the  end  of  every  bit  interval  by 
appropriately  staggering  the  hold  commands.  At  a  given  instant  of 
time  t  satisfying  mTb  <  t  <  {m+l)Tb,  stored  samples  of  the  integrator 
outputs  for  the  (m-l)-th  and  the  (ml-^h  bit  intervals  (e-g.,  exm, 

ex(m-l)*  eym*  ant*  ey(m-l)^  are  ava^a^^e  at  Ionr  of  the  six  S/H  outputs, 

so  that  the  dot  product  7tm-l)> =  exmex(m-l)  +  eymey(m-l)  can  be 

formed  in  the  multiplier  circuit.  Eight  analog  switching  circuits  are 

used  to  control  the  order  in  which  the  S/H  outputs  are  multiplied,  as 
*  ^ 

shown  in  Table  2  for  the  4>x  component;  the  arrangement  is  the  same 

A 

for  the  4>y  component. 

The  products  exrnex(m_i)  and  eymey^.i)  are  formed  by  two 
four-quadrant  multipliers.  At  the  multiplier  outputs,  the  message 
carrying  component  of  the  signal  is  a  binary  waveform  modulated  by  a 
sinusoid  having  a  frequency  equal  to  2AfQ.  Ideally,  the  sinusoidal 


TABLE  2 

The  cyclic  arrangement  of  the  analog  switching 
and  sample- hold  functions 


Signalling 

Interval 

Sample-Hold  Inputs 

Analog  Switches 

Sample-Hold 

Inputs  Multiplied 

#1 

#2 

#3 

#1 

#2 

#3 

#4 

(Multiplier 

Output) 

m 

HOLD 

HOLD 

SAMPLE 

OH 

OFF 

ON 

OFF 

#1  x  #2 

ro+1 

SAMPLE 

HOLD 

HOLD 

OFF 

ON 

OFF 

ON 

#2  x  #3 

bH-2 

HOLO 

SAMPLE 

HOLD 

ON 

OFF 

OFF 

ON 

#1  x  13 

m+3 

HOLD 

HOLD 

SAMPLE 

ON 

OFF 

ON 

OFr 

#1  x  #2 

m+4 

SAMPLE 

HOLD 

HOLD 

OFF 

ON 

OFF 

ON 

#2  x  #3 

m+5 

HOLD 

SAMPLE 

HOLD 

ON 

OFF 

OFF 

ON 

#1  x  #3 

modulation  on  the  output  of  one  multiplier  is  180*  out  of  phase  with 
respect  to  the  modulation  on  the  second  multiplier's  output.  The  two 
multiplier  outputs  are  added  using  an  (inverting)  operational  amplifier 
to  generate  Ideally,  the  sinusoidal  modulation  (at  the 

frequency  2AfQ)  at  the  amplifier  inputs  does  not  appear  at  the  output. 

The  signal  representing  Zm*Z^m_^  is  then,  applied  to  a  level  detector 
to  generate  the  detected  bit  stream.  If  the  input  to  the  level  detector 
is  greater  than  zero,  the  output  of  the  level  detector  is  a  positive  volt¬ 
age;  otherwise,  the  output  i3  zero. 

The  experimental  receiver  was  designed  to  give  satisfactory  per¬ 
formance  for  data  bit  rates  up  to  500  Kbps.  However,  at  high  bit  rates, 
demodulator  performance  is  degraded  by  bandlimiting,  slew  rate  limit¬ 
ing,  integrator  dump  time  limitations,  and  other  extraneous  receiver 
imperfections  (see  Appendix  C).  At  low  bit  rates,  the  time  required  to 
obtain  an  accurate  measurement  of  the  bit  error  probability  is  large. 
For  these  reasons,  the  experimental  data  to  be  presented  was  obtained 
with  the  rate  set  at  an  intermediate  value,  50  Kbps.  A  list  of  receiver 
parameters  for  f^  =  50  Kbps  is  given  in  Appendix  D. 

C.  Description  of  Peripheral  Test  Circuits 

A  block  diagram  of  the  circuits  used  to  bench  test  the  experi¬ 
mental  differential  detector  is  si: own  in  Figs.  13a  and  13b.  These 
circuits  were  resigned  to  perform  the  following  functions; 
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(1)  generation  of  the  additive  channel  noise  process,  the  bit  timing 
error  process,  the  data  bit  stream,  and  the  transmitted  waveform, 

(2)  measurement  of  the  bit  timing  error,  the  frequency  offset,  the 
signal  to  noise  ratio,  the  BEP  and  the  probability  of  two  consecutive 
errors,  and  (3)  computation  of  o^2  and  •  •  The  basic  operation  of  the 
peripheral  circuitry  is  described  in  this  section. 

The  bit  timing  jitter  and  the  channel  noise  processes  were  de¬ 
rived  from  separate  noise  sources  to  ensure  their  6tatist;cal  indepen¬ 
dence.  Each  source  produced  zero  mean  Gaussian  noise  in  a  band 
7  MHz  wide  centered  at  30  MHz.  The  spectra  of  the  two  noise  voltages 
was  flat  to  within  i 0.5  dB  over  the  frequency  range  30  MHz  ±1  MHz. 

A  hybrid  combiner  was  used  to  add  the  channel  noise,  n(t),  to  the 


Fig.  13a.  A  block  diagram  of  the  peripheral  test 
circuits  used  to  generate  the 
transmitted  waveform. 


K 

\ 
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Fig.  13b.  A  block  diagram  of  the  peripheral  test  circuits 
used  to  generate  receiver  imperfections  and 
to  measure  receiver  performance. 


transmitted  waveform  to  obtain  the  signal  Z'(t).  Since  2  MHz  »  50 
Kbps,  the  noise  voltage  approximated  a  white  Gaussian  process.  Due 
to  a  limitation  in  the  dynamic  range  of  the  mixer  stage  at  the  DD  input, 
the  signal  Z'{t)  was  bandlimited  to  Z  MHz(=Bj);  the  filter  output  was  the 
receiver  input  Z(t).  The  baseband  noise  used  to  generate  bit  timing 
error  was  obtained  by  down-converting  the  appropriate  30  MHz  noise 
signal  and  filtering  the  resulting  waveform  to  the  desired  bandwidth. 

The  data  bit  stream  consisted  of  a  maximum  length  127-bit  P-N 
sequence.  The  data  bit  stream  was  differentially  encoded  using  the 
circuit  diagrammed  in  Fig.  14  and  then  used  to  biphase  modulate  the 
30  MHz  carrier.  In  the  theoretical  analysis  of  the  effects  of  bit  timing 
error,  each  four  bit  sequence  was  assumed  to  occur  in  the  encoded  bit 
stream  with  equal  likelihood.  Using  the  127-bit  P-N  sequence  resulted 
in  an  approximation  to  this  assumption,  as  illustrated  by  Table  3.  The 
probability  of  occurrence  of  the  {0,  0,  0,  0}  sequence  is  equal  to  about 
0.551,  while  the  probability  of  occurrence  of  any  other  four  bit 
sequence  equaled  about  0.063.  Since  bit  timing  error  does  not  affect 
performance  when  {0,  0,  0,  0}  is  transmitted,  the  measured  bit  error 
probability  should  be  slightly  larger  than  the  value  determined  theore¬ 
tically.  The  effect  is  very  small,  however,  since  one  less  {0,0,  0,0} 
sequence  per  127  sequences  amounts  to  1  «.  than  one  percent. 
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Fig.  14.  A  block  diagram  of  the  differential 
encoder. 


TABLE  3 

The  probability  of  occurrence  of  each 
four  bit  sequence 


Sequence 

Probability 
of  Occurrence 

0000 

0.0551 

0001 

0.0630 

1000 

i; 

nn 

ii 

1110 

n 

0111 

ii 

0110 

ii 

0010 

ii 

0011 

ii 

1010 

ii 

ion 

ii 

1101 

ii 

1100 

ii 

0101 

ii 

0100 

ii 
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The  system  clock,  a  50  KHz  sinusoid,  is  used  to  generate  the 
absolute  system  clock,  C^,  the  phase  adjustable  clock,  Cp,  and  the 
phase  adjustable  timing  jitter  clock,  Cj.  The  absolute  system  clock, 
a  square  wave  derived  from  the  zero  crossings  of  the  50  KHz  system 
clock,  is  used  to  generate  an  error-free  data  bit  stream  for  compari¬ 
son  with  the  receiver  output.  The  phase  adjustable  system  clock  is 
used  to  generate  the  transmitted  bit  stream;  the  phase  of  Cp  is  adjust¬ 
ed  to  compensate  for  circuit  delays  with  respect  to  C^.  The  second 
adjustable  system  clock,  Cj,  is  used  to  time  the  digital  control  cir¬ 
cuits  in  the  DD;  its  phase  directly  determines  the  bit  timing  offset  *n. 

The  bit  timing  error  is  generated  by  comparing  a  baseband  noise 
voltage  with  a  sawtooth  waveform.  The  period  of  the  sawtooth  equals 
the  data  bit  length,  T5,  as  shown  in  Fig.  15,  and  its  phase  is  deter¬ 
mined  by  Cj.  The  noise  voltage,  Vj,  is  a  narrowband,  zero-mean 
Gaussian  process  with  a  bandwidth  equal  to  Z  KHz.  The  sawtooth  and 
noise  voltages  are  compared  using  an  analog  comparator.  The  output 
of  the  comparator  is  a  series  of  pulses  whose  leading  edge  remains  in 
phase  with  the  sawtooth  and  whose  trailing  edge  is  delayed  t'  seconds 
with  respect  to  the  leading  edge  (see  Fig.  15): 


t'  = 


ZVA 


rb  +  to 


(95) 


Fig.  15.  Generation  of  bit  timing  jitter. 


If  Vt  is  assumed  constant  during  each  period  of  the  sawtooth  waveform, 


then  it  can  be  shown  that  t'  is  a  sample  function  from  a  Gaussian  pro¬ 


cess  with  mean  and  variance  given  by 


(96) 


E(t')  =  t 


and 


(97)  <rjt  = 

t  4  VA2 
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The  timing  of  the  integrate/ dump  commands  are  derived  from  the 
timing  of  the  trailing  edge  of  each  pulse  at  the  comparator  output;  thus 
<rt*  is  the  variance  of  the  bit  timing  error  process.  The  variance  <r€* 
is  varied  by  adjusting  the  variance  of  Vj,  and  the  bit  timing  offset,  «, 
is  varied  by  adjusting  the  phase  of  the  sawtooth  waveform  with  respect 
to  the  timing  of  the  received  signal  Z(t),  i.e.»  the  phase  of  Cj  with 
respect  to  Cp. 

Since  V*  is  only  approximately  constant  during  each  decision  in¬ 
terval,  the  length  of  the  integration  interval  varies  slightly  with  time. 
To  see  that  this  deviation  is  negligible,  consider  the  following  example. 
Suppose  that  the  normalized  standard  deviation,  <r<n(*  «r(/T^),  is  equal 
to  0.10.  Then  most  values  |«n  -  tn|  will  be  within  3ir<n.  If  a  sinusoidal 
time  variation  of  <  at  a  frequency  of  2  KHz  and  an  amplitude  of  3 or,  is 
assumed,  then  the  maximum  variation  in  the  length  of  an  integration 
interval  will  be  equal  to  about  (0.01)Tb»  assuming  a  50  Kbps  data  bit 
rate.  This  results  in  a  ±  1%  variation  in  the  effective  signal  power  to 
noise  power  ratio  at  the  integrator  output,  a  very  small  deviation. 

The  average  variations  in  the  integration  interval  w'll  be  much  smaller 
than  +  1%.  Since  the  maximum  normalized  standard  deviation  consid¬ 
ered  in  this  report  will  be  0,10,  the  effects  of  variations  in  the  inte¬ 
gration  interval  will  be  negligible. 
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A  time  interval  meter  (TIM)  is  used  in  conjunction  with  a  digital 
computer  to  determine  the  statistics  of  the  bit  timing  error.  The 
sta^c  and  stop  commands  required  to  control  the  TIM  are  generated 
using  circuits  having  the  block  diagrams  shown  in  Fig.  16.  Start 
pulses  are  generated  from  the  absolute  system  clock  at  a  rate  of  six 
pulses  per  second.  The  stop  pulses  are  generated  by  digitally  dividing 
the  frequency  of  Cj  by  two.  The  TIM  measures  the  time  interval  be¬ 
tween  the  start  and  stop  pulses.  Since  the  TIM  can  read  only  positive 
time  intervals,  the  stop  pulse  is  delayed  20  ^sec  with  respect  to  the 
first  pulse.  The  bit  timing  error  is  thus  related  to  the  TIM  reading 
(A  t)  by 

(98)  <  =  At  -  20  x  10'6  +  D  . 


Fig.  16,  A  block  diagram  of  the  bit  timing 
error  sampler. 
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The  constant  D  takes  into  account  delays  introduced  by  the  digital  and 
analog  circuits.  Its  value  is  determined  by  synchronizing  the  inte¬ 
grator  dump  command  with  the  bit  transitions  at  the  integrator  input, 
a  condition  that  corresponds'to  a  zero  bit  timing  error.  It  was  found 
that  this  point  could  be  determined  within  an  accuracy  of  +  20  nsec. 

The  TIM  is  accurate  within ,+  10  nsec.  The  mean  («),  variance  (O'2), 
and  the  histogram  of  the  bit  timing  error  process  were  estimated  from 
3000  samples  of  the  bit  timing  error.  Histograms  for  various  bit 
timing  jitters  showed  that  the  bit  timing  error  process  closely  approx¬ 
imated  a  Gaussian  process. 

The  signal  to  noise  ratio  was  measured  using  the  30  MHz  IF 
stage  of  a  receiver.  The  IF  bandwidth  is  nominally  equal  to  5  KHz  and 
is  stable  to  within  +0.5%.  The  power  of  the  message  carrying  com¬ 
ponent  of  the  received  signal*  Pr.  was  measured  by  removing  the  bi¬ 
phase  modulation  and  the  additive  channel  noise.  Similarly,  the  noise 
power  was  measured  (in  the  5  KHz  bandwidth)  by  reducing  the  trans¬ 
mitted  power  to  zero.  The  bit  energy  to  (single-sided)  noise  density 
ratio,  E^/Nq,  was  then  calculated  from  these  measurements. 

Two  error  probabilities,  Pg  and  PEjy  were  estimated  by  aver¬ 
aging  each  error  count  over  a  large  number  of  bit  intervals.  Each 
measurement  was  entered  as  a  data  point  if:  1)  the  error  count  ex¬ 
ceeded  100,  2)  the  number  of  observed  bit  intervals  exceeded  3x  10*, 
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and  3)  the  point  was  repeatable  to  within  ±1%  with  respect  to  the  bit 
energy  to  noise  density  ratio. 


CHAPTER  V 

EXPERIMENTAL  RESULTS 


A.  Best  Case  Performance 

The  bit  error  probability  of  the  experimental  receiver  was  mea¬ 
sured  for  the  Afc  »  0,  <  =  0,  6  =  350  nsec  case  to  determine  the  effects 
of  extraneous  receiver  imperfections  not  included  in  the  mathematical 
model*  These  effects  include  bandlimiting,  circuit  non-linearities, 
d.c.  offsets,  non-zero  aperture  and  acquisition  times  in  the  sample/ 
hold  circuits,  error  in  matching  the  gain  in  each  orthogonal  channel, 
error  in  LO  phasing,  and  measurement  error.  If  all  receiver  imper¬ 
fections  could  be  eliminated,  the  bit  error  probability  versus  E|,/N0 
would  be  given  by  Eq.  (23),  The  measured  results,  compared  to  this 
ideal  result,  are  shown  in  Fig,  17,  Each  data  point  was  within  0,17  dB 
of  ideal.  From  Eq.  (24),  it  can  be  shown  that  about  0.075  dB  of  loss 
is  due  to  a  finite  integrator  dump  time.  The  remaining  0.095  dB  loss 
is  primarily  due  to  the  effects  of  the  finite  bandwidth  Bj  and  bandlimit¬ 
ing  in  the  integrate/dump  and  sample/hold  circuits  (see  Appendix  C). 
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Fig.  17.  The  bit  error  probability  versus  the  bit  energy 

to  noise  density  ratio  under  best  case  conditions; 
«*0,  (rf*0  and  Afo*0.  Both  analytical  and  experi¬ 
mental  results  are  given.  Analytical  and  experi¬ 
mental  results  for  the  probability  of  two  conse¬ 
cutive  errors  are  also  given. 
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B.  Double  Error  Measurements 

The  probability  of  error  given  that  the  demodulator  erred  in  the 
previous  bit  interval,  Pr(em^em-1^*  was  determined  experimentally 
as  a  function  of  E^/No  for  «n  »  0,  r*n  s  °»  and  Af0  *  0.  The  results, 
given  in  Fig.  18,  compared  extremely  well  with  the  theoretical  result 
for  E5/ N0  in  the  range  4  to  9  (numeric).  The  deviation  from  the 


'  (  NUMERIC) 


Fig.  18.  The  probability  of  error  given  that  an  error 
was  detected  in  the  previous  bit  interval 
versus  the  bit  energy  to  noise  density  ratio 
for  «=0,  0^=0  and  Afo=0;  both  analytical  and 
experimental  results  are  given. 
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theoretical  curve  is  somewhat  larger  for  the  data  point  corresponding 
to  Efc>/N0  =  11.2  because  the  ratio  PEq/Pe  *s  extremely  sensitive  to 
receiver  imperfections  such  as  intersymbol  interference  and  band- 
limiting;  that  is,  at  high  Eb/N0,  the  ratio  of  the  theoretical  and  mea¬ 
sured  BEP  becomes  greater  even  though  the  performance  difference  in 
terms  of  E}j/N0  remains  approximately  the  same  (see  Fig.  17). 

C.  The  Effects  of  Bit  Timing  Error 

Curves  showing  the  bit  error  probability  as  a  function  of  the  bit 
energy  to  noise  density  ratio  for  selected  combinations  of  bit  timing 
offset  and  bit  timing  jitter  are  given  in  Figs.  19-22;  both  analytical  and 
experimental  results  are  shown.  Again,  the  analytical  and  experi¬ 
mental  results  are  in  close  agreement. 

It  is  recalled  that  the  results  of  the  theoretical  analysis  showed 
that  the  bit  error  probability  was  invariant  to  the  sign  of  the  bit  timing 
error.  Referring  to  Fig.  19,  it  is  seen  that  the  measured  bit  error 
probability  depended  on  the  sign  of  the  bit  timing  error.  This  result 
was  due  to  bandlimiting  in  the  integrate/dump  and  sample/hold  cir¬ 
cuits,  which  had  the  effect  of  introducing  a  171  nsec  delay  at  the 
sample/hold  output  in  response  to  a  step  applied  at  the  integrator  in¬ 
put  (see  Appendix  C).  The  effect  o.f  the  delay  is  illustrated  in  Fig.  23. 
When  the  bit  timing  error  is  positive,  the  overall  effect  of  the  delay  is 
to  reduce  the  signal  loss  that  an  ideal  integrator  and  sampler  would 
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experience  with  bit  timing  error.  For  a  negative  bit  timing  error,  the 
delay  causes  the  signal. loss  to  be  greater.  Consequently,  the  BEP  was 
smaller  for  positive  bit  timing  offsets  than  for  negative  bit  timing  off- 
sets. 

Results  given  in  Figs.  20-22  for  cases  where  both  the  mean  and 
the  variance  are  non-zero  apply  when  the  timing  offset  is  negative; 
thus,  they  show  the  detector's  performance  under  worst  case  condi¬ 
tions.  The  BEP  for  7n  positive  was  slightly  lower.  Note  that  for  the 
7n  =  0,  ir<n  =  0.10  and  the  7n  =  °*  r«n  =  cases,  the  measured  BEP 
represents  the  average  effect  of  both  positive  and  negative  offsets. 

It  should  also  be  noted  that,  when  the  bit  timing  error  is  non¬ 
zero,  receiver  performance  may  be  improved  by  increasing  the  dump 
time  6.  For  example,  consider  the  effect  of  dump  time  when  «n  = 

+  0.08.  If  6  is  set  equal  to  0.16  Tb»  then  the  bit  timing  offsets  having 

magnitudes  less  than  0.08Tb  do  not  affect  the  message  carrying  com- 

Eb 

ponent  of  the  integrator  outputs.  Correspondingly,  for  *  10  dB,  the 


bit  error  probability  is  given  by 


Now,  the  bit  error  probability  for  the  otherwise  ideal  differential  de¬ 
tector,  with  7n  =  0.08,  equals  2.08  X  10"'“--  twice  as  large  as  the 
result  in  Eq.  (99).  In  certain  instances,  then,  the  imperfectly-timed 
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BIT  ERROR  PROBABILITY 


Fig.  19. 


The  bit  error 
noise  density 
«n=±0 .08  and 
experimental 


probability  versus  the  bit  energy  to 
latio  for  the  imperfectly-timed  DD; 
~n=±0.16.  Both  analytical  and 
results  are  given. 
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BIT  ERROR  PROBABILITY 


Fig.  20.  The  bit  error  probability  versus  the  bit  energy  to 
noise  density  ratio  tor  the  imperfectly-timed  DD; 
«n=0  and  <rsn=0.06.  Both  analytical  and  experi¬ 
mental  results  are  given.  Nomi:  al  experimental 
results  for  the  probability  of  two  consecutive 
rrror'i  are  .ilsn  given. 
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BIT  ERROR  PROBABILITY 


Fig.  21.  The  bit  error  probability  versus  the  bit  energy  to 
noise  density  ratio  for  the  imperfectly-timed  DD; 
Tn=0  and  <r<n=0.10.  Both  analytical  and  experi¬ 
mental  results  are  given.  Nominal  experimental 
results  for  the  probability  of  two  consecutive 
errors  are  also  given. 
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BIT  ERROR  PROBABILITY 


«n=0.06  and  7n=- 0,08.  Both  analytical  and  exper¬ 
imental  results  are  given.  Nominal  experimental 
results  for  the  probability  of  two  consecutive 
pvn'i's  avo  alpo  given. 
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DD  will  have  a  lower  bit  error  probability  if  the  dump  time  i3  set  to 
some  non- zero  value. 


D.  The  Effects  of  Frequency  Offset 

Curves  comparing  the  measured  and  analytical  BEP  as  a  function 


(■ 


of  Eb/N0  for  selected  values  of  the  normalized  frequency  offset;  A  f„ 

A  <*>oTb  \ 

=  AfQTb  = — - - j  ,  are  given  in  Fig.  24.  Some  experimental  and 

analytical  results  are  also  compared,  in  terms  of  the  ratio  p/ pQ  (see 
p.  38)  as  a  function  of  Eb/N0,  in  Fig.  25.  All  experimental  results 
were  found  to  closely  agree  with  analytical  results  for  both  positive 
and  negative  offsets.  In  addition,  spot  measurement  results  of  the 
BEP  for  frequency  offsets  in  the  range  ir/2  ^  |  Aw0Tb  |  <  3u7 2  closely 
agreed  with  the  bit  error  probability  Pg  given  by 


(100) 


Pe(Au0,  Li)  =  u(|Au0Tb|  -it / 2 )  -  sgn(  | Aw0Tb  |  -ir/2)  PE 


where 

(101) 


x  <  0 
x  >  0 


(102)  sgn(x)  =  2 


-1  x  <  0 
+  1  x  >  0 


$ 

Some  of  the  results  shown  in  Fig.  25  a*-e  discussed  in  the  following 
section. 
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ig.  24.  The  bit  error  probability  versus  the  bit  energy 
to  noise  density  ratio  for  selected  values  of 
frequency  offset,  both  analytical  and  experi¬ 
mental  results  are  given. 
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J 


t 


1 

and  Pg  is  given  by  Eq.  (94),  even  though  the  assumption  Aw0Tb  <  ir/ Z  \ 

made  in  the  derivation  of  Eq.  (94)  is  no  longer  valid.  This  is  due  to  j 

j 

the  fact  that  the  double  frequency  components  in  4x(t)  and  £y(t)  are  j 

highly  attenuated  by  bandlimiting  in  the  integrate /dump  and  the  sample/ 
hold  circuits,  a  factor  not  taken  into  account  in  the  derivation, 

£,  Additional  Experimental  Results 

The  combined  effects  of  frequency  offset  and  bit  timing  error  on 
the  bit  error  probab  ility  and  the  probability  of  two  consecutive  errors 
are  extremely  difficult  to  investigate  analytically;  for  this  reason,  an 
experimental  investigation  was  performed.  Since  the  experimental 
results  previously  given  agreed  quite  well  with  available  analytical 
results,  the  results  to  be  given  should  accurately  reflect  receiver  per¬ 
formance.  Results  were  obtained  showing  1)  the  combined  effects  of 
bit  timing  error  and  frequency  offset  on  the  BEP,  2)  the  effect  of  bit 
timing  error  on  the  probability  of  two  consecutive  errors,  and  3)  the 
effect  of  frequency  uncertainty  on  the  probability  of  two  consecutive 
errors . 

The  combined  effects  of  bit  timing  error  and  frequency  offset  on 
the  BEP  as  a  function  of  Eb/No  was  determined  for  selected  values  of 
«n;  the  results  are  given  in  Figs.  26  and  27.  The  performan:e  degra¬ 
dation  due  to  frequency  offset  for  the  perfectly  and  imperfectly  timed 
detectors  can  be  compared  by  redefining  the  ratio  p/p0  (see  p.  38)  to 
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Fig.  26.  The  bit  error  probability  versus  the  bit  energy  to 
noise  density  ratio  of  the  imperfectly-timed  DD 
for  selected  values  of  the  frequency  offset;  nominal 
experimental  results  are  given  for  the  ~in  =  -  ^.U8 
case . 
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Fig.  27.  The  bit  error  probability  versus  tin  bit  energy  to 
noise  density  ratio  of  the  imperfectly-timed  DC 
for  selected  values  of  the  frequency  offset; 
nominal  experimental  results  are  given  for 
the  7n  =  -  0 . 16  case . 
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include  the  effects  of  bit  timing  error.  In  Fig.  25,  measurements  of 
the  ratio  are  compared  to  theoretical  results  for  the  Tn=  -  0.16  . 
and  the  «n=  -  0.08  cases  with  Pjr=  10"4.  In  this  figure,  ji0  and  p,  have 
been  defined  as  the  bit  energy  to  noise  density  ratios  required  to  main¬ 
tain  a  given  BEP  at  the  output  of  an  imperfectly-timed  DD  for  a  zero 
and  non-zero  frequency  offset,  respectively.  Note  that  the  normaliza¬ 
tion  factor  is  not  the  same  for  different  bit  timing  offsets,  even 
though  the  BEP  is  the  same.  The  results  show  that  the  performance 
degradation  caused  by  a  given  frequency  offset  decreases  only  slightly 
even  when  the  bit  timing  offset  is  large,  e.g.,  7n=-  0.16.  This  indi¬ 
cates  that  the  ratio  |jl/ p0  for  the  imperfectly-timed  demodulator  can  be 
estimated  from  the  theoretical  curves  given  in  Fig.  11. 

The  probability  of  two  consecutive  errors  was  measured  for 
selected  combinations  of  bit  timing  jitter  and  bit  timing  offset;  results 
are  given  in  Fig.  28  in  the  form  of  the  ratio  Pr(emlem-l)  = 
versus  E^/Nq.  The  results  indicate  that  increases  in  uncertainty  in 
the  bit  arrival  time  leads  to  increases  in  the  consecutive  error  ratio.. 
This  is  a  reasonable  result,  since  bit  timing  error  introduces  an 
effective  signal  to  noise  ratio  loss  at  the  integrator  outputs.  Again, 
it  was  found  that  Pr(em  |  em_  j )  decreased  as  the  bit  energy  to  noise 
density  ratio  was  increased,  except  for  the  7n  =- 0*  16  case.  For 
most  cases  of  practical  interest,  Pr(em  |  em_  j)  will  be  less  than  25%. 
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Fig.  28.  The  probability  of  emr  given  that  an  error 
was  detected  in  the  previous  bit  interval 
versus  the  bit  energy  to  noise  density  ratio 
for  the  imperfectly-timed  DD.  nominal 
experimental  results  are  given. 
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The  probability  of  an  error  given  that  a  detection  error  occurred 
in  the  previous  bit  interval  as  a  function  of  frequency  offset  was  mea¬ 
sured  for  selected  values  of  the  bit  energy  to  noise  density  ratio; 
results  are  given  in  Fig.  29.  It  was  found  that  the  ratio  Pr(emlem-l) 
decreases  rapidly  as  the  frequency  offset  is  increased  from  zero  uuHl 
a  minimum  is  reached  at  Afn  -  0.12  =  0.75  rad).  The  mea¬ 

sured  BEP  versus  Afn  is  also  shown  in  Fig.  29.  Note  that,  for 
l^nl  >  0*09,  the  ratio  Pr(em  |  em_i )  is  less  than  the  bit  error  pro¬ 
bability;  that  is,  for  normalized  frequency  offsets  greater  than  0.09, 
an  error  is  less  likely  to  occur  if  a  detection  error  occurred  in  the 
previous  bit  interval.  It  was  also  found  that  the  probability  of  two  con¬ 
secutive  errors  (PEj))  decreased  slightly  with  increasing  frequency 
offset  to  a  point,  as  ’ndicated  in  Fig*  30.  For  the  values  of  E^/N0 
tested,  Pejj  reached  a  minimum  value  at  a  normalized  frequency  off¬ 


set  of  about  0.10 


Fig.  29.  The  probability  of  error  given  that  an  error  wa 
cl  tecied  in  the  previous  bit  interval  versus  fre 
quency  offset  for  selected  values  of  the  bit 
energy  to  noise  density  ratio;'  nominal  experi¬ 
mental  results  a*e  given.  Analytical  results 
of  the  bit  error  probability  fc'  each  case  are 
also  given. 
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Fig.  30.  The  probability  of  two  consecutive  errors 

versus  frequency  offset  for  selected  values 
of  the  bit  energy  to  noise  density  ratio; 
nominal  experimental  results  are  given. 
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CHAPTER  VI 

SUMMARY  AND  CONCLUSIONS 


An  extensive  analytical  and  experimental  investigation  of  DPSK 
and  differential  detection  has  been  presented.  The  results  provide  a 
basis  for  the  design  and  implementation  of  a  practical  differential 
detector  in  which  changing  system  geometries  are  involved,  ;. g.,  in  a 
satellite  communications  link.  Available  analytical  results  relating  to 
the  performance  of  the  ideal  and  imperfectly-timed  differential  detec¬ 
tor  have  been  summarized  and  an  expression  for  the  bit  error  probabil¬ 
ity  in  cases  where  the  estimate  of  the  carrier  frequency  is  in  error  h;.s 
been  derived  and  numerical  results  documented.  To  check  the 
accuracy  of  analytical  results,  an  experimental  differential  detector 
was  implemented.  The  details  of  the  detector's  construction  have  been 
described  and  its  performance  documented.  The  experimental  results 
were  found  i  closely  agree  with  available  analytical  results  and  with 
the  result  derived  relating  to  the  effects  of  frequency  uncertainty.  It 
is  concluded  that  the  analytical  results  presented  in  this  report  can  be 
applied  to  practical  DPSK  systems  with  confidence  and  that  the  experi¬ 
mental  results  obtained  for  cases  in  which  no  analytical  results  are 
available  accurately  reflect  demodulator  performance. 


APPENDIX  A 

A  COMPARISON  OF  IDEAL  COHERENT  AND 
IDEAL  DIFFERENTIAL  DETECTORS 


If  the  additive  channel  noise  is  a  zoro  mean  hite  Gaussian  pro¬ 
cess  and  antipodal  modulation  is  employed,  then  the  bit  error  pro¬ 
babilities  for  ideal  PSK,  ideal  DC-PSK,  and  ideal  DPSK  (using  differ¬ 
ential  detection)  are  given  by 
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for  DPSK[6,  7,  9] 


Curves  showing  the  additional  signal  energy  required  for  ideal  differ¬ 
ential  detection  to  maintain  a  BEP  equal  to  that  of  an  ideal  coherent 
detector  and  for  an  ideal  differentially-coded  coherent  detector,  as  a 
function  of  the  bit  energy  to  noise  density  ratio  at  the  input  to  the  co- 
herent  detector,  is  given  in  Fig.  31.  At _ _  =  10  dB,  the  ideal  coher- 

N0 

ent  detector  performs  only  0.70  better  than  the  ideal  DD.  This  is  sign 
lficant  in  that  the  DD  does  not  require  a  phase  reference. 
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Fig.  31.  A  comparison  of  ideal  differential  detection,  'deal 
differentially-coded  coherent  detection,  and  ideal 
coherent  detection. 
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APPENDIX  B 
CALCULATIONS 


The  following  expression  for  Pg  is  to  be  simplified  to  facilitate 
numerical  evaluation: 
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The  integral  E2  will  first  be  evaluated.  It  can  easily  be  shown  that 
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Therefore,  E2  is  invariant  to  the  sign  of  the  frequency  offset.  Also, 
it  is  easy  to  show  that 
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(112) 
(cont. ) 
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cos  AwcTb 


[-  t]  *> 

i  LjWjcI  tosx 

L  No  . 

[‘  r] du 


.  exp  [-  Jb!i !!»!* 

L  N0  J 


cos  x  dx 


By  using  the  relation[l] 


(113) 


W2N^  ‘'O  ^0 


^ir/Z  „B  —  cos  i|/i  f  r 
V  \  4N°  exp<-  - 

J0  J0  L  L 


Z  No  JJ  dljj. 


;  m  r’ 
'  2  '  it  I 


-a2b2/nd 

exP  A2TosTeTBrsh7e' 


and  by  ietting  A  =  Lj  and  B  *  Lj  cosAw0Tb  =  l2  ,  the  first  term  in 
Eq.  (112)  reduces  to 


(114) 


U  -  cos  x  cos  Aw0Th 

2l,  r2r  Jno 

rrzr  4)  4 


W3n^ 


(  _/u2  +  LiZ  sin2x^ 
L\2  N0  y_ 


•  cosxdx 


1  1  rr/z 

=  2'«J0  exp;' 

1  >  fW2  r 

expb 


V  Li2  cos2Au>0Tv 


eXP  Li2cos29+  Li2 coszAwQTb  sin2 9 


L,2L  2 


)0  exp  L“  L12cos're+  L,2sm2oJ  d° 


Consequently, 


that  Ei ,  given  by 


(117) 


it/ 2  , 


El=ikf‘A  f 

•s/ttNq  Jq  ^0 


Lj  - cos  i|> 

'  jNo 


simplifies  to 


(118)  E,  =  exp  - 


Thus, 


(119) 


!xp["  7"]du 

|  _  Lx2sin2<]itl 

L  nq  J 


cos  4»i 


■if 


Li*L;2/N0 

Lj2  cos2 9+  L22sin^0 


Lx  - cos(4i+Aw  Tb) 

_ur/2fK  «Pr4i 

"^o  -W2  L,[Z  coS  + cos A%Tb  J 

JNo 


Lx2  sin2 4^1 
f  No 


cos  4*  dijj< 


APPENDIX  C 

THE  INTEGRATE/ DUMP  AND  SAMPLE/ HOLD 

CIRCUITS 


A .  Tne  Integrate/ Dump  Circuit 

A  schematic  diagram  of  the  integrate/dump  circuit  employed  in 
the  DD  is  given  in  Fig.  32.  The  circuit  was  designed  to  provide  accu¬ 
rate  integration  for  small  valuec  of  time  and  to  have  a  very  short  dump 
time  capability.  To  evaluate  actual  perfoxTnance,  the  integrator  will 
be  analyzed  for  the  short  integration  interval  case. 


OP 

AMP 


15V.  INTEGRATEr~2N4859 
OV.  DUMP  ’-PET 


Fig.  32.  A  schematic  diagram  of  the  integrate/ dump 
circuit  employed  in  the  DD. 


If  She  operational  amplifier  in  Fig  3?,  <;  assumed  to  be  ideal, 


e.g.,  if  infinite  gain,  bandwidth,  and  input  impedance  is  assumed, 
then  the  output  voltage  VQ(t)  is  given  by 

T 

(120)  yo(t)  -  A  f  v.  dt  . 

In  practice,  the  open  loop  gain,  bandwidt.  ,  and  input  impedance  are 
finite.  For  example,  the  transfer  function  of  most  operational  ampli¬ 
fiers  has  a  6  dB/octa”e  rolloff  at  high  frequencies,  as  illustrated  in 
Fig.  33.  On  assuming  that  the  input  impedance  is  infinite,  that  the 
input  bias  current  and  voltage  offset  are  negligible,  and  that  the  ampli¬ 
fier  has  the  transfer  function  given  in  Fig.  33,  it  can  bo  shown  that 

V 

the  transfer  function _ 2  (S)  is  given  by 


(121) 


Ro 

AR  +  1 


A0R 


+  — +  ^+~p  +  JR£_S* 

irf-p  AQ  Z  J  Zirf-p 


whire  Aq  represents  the  open  loop  d.c.  gam,  fj  the  small  signal 
unity  gain  bandwidth,  and  RQ  the  open  loop  output  resistance.  If,  in 
addition,  it  is  assumed  that 


(122)  - <  <  <  1 
AtR 


8b 


(123) 


irRCfy  >  >  >  1 


Fig.  33.  Bode  diagram  comparing  the  transfer  function 
an  ideal  integrator  to  that  of  a  practical 
integrator . 


and 


(124)  Aq>  >  >  1 


then  Eq.  (121)  can  be  approximated  by[M] 


H  7 


(125) 


1 


A  Bode  diagram  of  the  above  expression,  given  in  Fig.  33,  shows  that 

the  practical  integrator  departs  from  an  ideal  integrator  for  large 

values  of  time  and  for  small  values  of  time. 

For  small  values  of  time,  TTR0Cf'p  will  normally  be  much  less 

than  the  open  loop  gain  Aq,  but  the  relationship  given  in  Eq.  (123)  may 

not  hold  if  _fL  approaches  (=  Z-n-f-r).  Consequently,  for  a  short  inte- 
RC 

V 

gration  period,  _ °{S)  is  approximated  by 

Vi 


(126) 


The  output  voltage  in  response  to  an  input  step  of  amplitude  E  is  given 
by 

(127)  VQ(t)  =  EF  {t  +  r0[  exp(-  V)  -  1  J  }  u(t) 


where 


(128) 


and 

(129) 


Zltiry 

F  =  - ± - 

2  +  rrRCf-p 


RC 

2  +  TtRCfT 
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Equation  (127)  shows  that  the  finite  bandwidth  fq-  introduces  a  time 
lag  equal  to  tq  seconds  at  the  integrator's  output,  as  illustrated  in 
Fig.  34.  Also,  if  the  constant  2(RC)  1  is  of  the  same  order  of  magni¬ 
tude  as  then  the  gain  of  the  integrator  will  depend  on  f  this  is 
usually  undesirable  because  of  temperature  instaoiiicies  or  parameter 
variations  among  op  amps. 

Substituting  circuit  parameters  of  the  elements  used  in  the  DD 
integrate/ dump  circuit  at  the  50  Kbps  rate,  given  by 
R  =  2KS2 


C  =  0 . 00  >  jif 


To 


Fig.  34.  The  step  response  of  an  integrator' 
lor  small  values  of  time. 


8l; 


and 


fT  =  10.  MHz 

into  Eqs.  (1.28)  and  (129),  it  is  found  that 

(130)  tq  =  31  nsec 

(131)  F  =  5  X  105 

and 

rTb  2  X  10‘5 

(132)  G  5  F  '  dt  =  5  x  105  \  dt  =  10 

C  J0 

At  a  500  Kbps  rate,  a  R,  G,  and  fj  held  constant,  td,  C,  and  F 
are  found  to  be 

(133)  tq  =  27  nsec 

% 

(134)  C  ^  168  pf 


and 

(135)  F  =  5  X  106 


The  gain  sensitivity  to  variations  in  fT  at  the  500  Kbps  rate, 


(136) 


1  _dF " 

FdfTL  . 


T=  10  MHz  fT(Z  +ttRC^t)J  fT  _  10MHz 


=  1 .58%/ MHz 


is,  for  most  applications,  acceptable. 
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The  minimum  dump  time  for  the  circuit  of  Fig.  3 2  is  deter¬ 
mined  by  1)  the  turn-on  time,  t^(0n)*  1  FET  switch,  2)  the  on 

resistance,  Ron,  of  the  FET  switch,  3)  the  value  of  the  integrating 
capacitor  (C),  4)  the  maximum  accept  Able  intersymbol  interference, 
and  5)  the  slew  rate,  Sr,  of  the  operational  amplifier.  At  the  50  Kbps 
rate,  the  circuit  parameters  were  given  by 

*on=25fi 
C  =  0.002  jif 

Sr  =  70  V/psec 
ld(on)  =  6  nsec 

The  dump  interval  was  set  so  that  the  integrator's  output  voltage  at 
the  end  of  the  dump  interval  equaled  0.002  times  the  output  voltage  ,it 
the  beginning  of  the  dump  interval;  correspondingly,  at  least  6.2  time, 
constants  must  be  spanned  by  the  dump  interval; 

(137)  ^MIN  **  ^ ^  (Ron  ^d(on)  ~  316  nsec  . 

Since  the  worst  case  rise  time  of  the  op  amp  is  about  175  nsec,  slew 
rate  limiting  does  not  need  to  be  considered  in  this  case.  However, 
at  the  500  Kbps  rate  , 

*Here,  worst  case  rise  time  is  defined  as  the  time  required  for  the 
op  amp  to  settle  to  within  0.002  times  the  full  output  voltage  after 
the  FET  has  turned  on  minus  the  delay  rQ. 


VI 


(138)  6.2{RonC)  -  26  nsec. 

Consequently,  the  minimum  dump  time  is  limited  by  the  slew  rate  to 
175  nsec  in  th’s  case.  At  either  bit  rate,  the  switching  time  of  the 
J-FET  is  relatively  small. 

The  dump  time  and  the  response  for  small  values  of  time  can  be 
improved  by  using  high  speed  op  amps.  For  example,  differential  op 
amps  with  slew  rates  of  1000  V/psec  and  f>p  =  100  MHz  are  presently 
available;  thus,  the  circuit  given  in  Fig.  32  can  provide  acceptable 
performance  for  bit  rates  of  5  Mbps  and  above  if  appropriate  devices 
are  utilized. 

B.  The  Sample/ Hold  Circuit 

A  schematic  diagram  of  the  sample/hold  circuit  is  given  in 
Fig.  35.  It  can  be  shown,  in  a  manner  similar  to  the  analysis  of  the 
integrate/ dump  circuit,  that  the  delay  t0'  of  the  output  voltage  in 
response  to  a  ramp  input  is  given  by 

(139)  tq'  =  120  nsec  . 

The  aperture  time,  defined  as  the  sample  to  hold  switching  delay,  was 
equal  to  about  20  nsec.  Therefore,  *he  to*"al  effective  delay  at  the  out¬ 
put  of  the  sample/ hold  circuit,  in  response  to  a  step  at  the  integrator's 
input,  is  gwen  by 

(140)  D0  =  x0  +  t  0  '  +  20  nsec  =  171  nsec 
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Fig.  35.  A  schematic  diagram  of  the  sample/hold 
circuit  employed  in  the  DD. 


APPENDIX  D 

A  LIST  OF  RECEIVER  PARAMETERS 


IF  BANDWIDTH 
DATA  BIT  RATE 

LO  PHASE  QUADRATURE  ERROR 
DELAY  AT  S/H  OUTPUT  IN  RESPONSE  TO 
A  STEP  APPLIED  AT  THE 
INTEGRATOR'S  INPUT 
S/H  APERATURE  TIME 
TOTAL  MULTIPLIER  ERROR 
MULTIPLIER  BANDWIDTH 
DATA  BIT  RATE  -  DESIGN  MAXIMUM 
ERRO.’'  IN  MATCHING  GAIN  IN  EACH 
ORTHOGONAL  CHANNEL 


2  MHz 
50  Kbps 
±  2° 


171  nsec 
c.0  nsec 

±  1%  FULL  SCALE 
3  MHz 
500  Kbps 


±  1% 
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